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BCS theory describes the formation of Cooper pairs and 
their instant "Bose condensation" into a superconducting 
state. Helium atoms are preformed bosons and, in addi- 
tion to their condensed superfluid state, can also form a 
quantum solid, lacking phase-coherence. Here we show 
that the fate of Cooper pairs can be more varied than the 
BCS or helium paradigms. In copper-oxide d-wave super- 
conductors (dSC) Cooper pairs are non-local objects, with 
both center-of-mass and relative motions. As doping de- 
creases, the center-of mass fluctuations force a correlated 
dSC into a state with enhanced diamagnetism and ro- 
bust but short-ranged superconducting order. At extreme 
underdoping, the relative fluctuations take over and two 
pseudogaps - "small" (charge) and "large" (spin) - emerge 
naturally from the theory, as Cooper pairs "disintegrate" 
and charge "detaches" from spin-singlet bonds. The en- 
suing ground state(s) are governed by antiferromagnetic 
rather than by superconducting correlations. The theory 
is used to account for recent experiments and to draw gen- 
eral conclusions about the phase diagram. 

Recent experiments [1] have narrowed the field of con- 
tenders for theoretical description of high-T c cuprates while 
simultaneously promoting the nature of the pseudogap state 
to the key conceptual issue. Two basic ideas are in play 
Q: either the pseudogap is a quantum disordered d-wave su- 
perconductor, or an entirely different form of a "competing" 
order, originating from the particle-hole (diagonal) channel 

EMU). 

In this paper, I show that the pseudogap in a correlated lat- 
tice dSC is both. Several recent experiments offer important 
clues in this respect: the observation of enhanced quantum 
diamagnetism in the pseudogap state of LSCO |6], alongside 
the giant Nernst effect |6], points to an intimate relation be- 
tween the pseudogap and a quantum disordered d-wave su- 
perconductor (dSC) |7] - there is hardly another known mi- 
croscopic mechanism which can deliver diamagnetism of this 
magnitude. Furthermore, the angle-resolved photoemission 
(ARPES) in underdoped LBCO |8] reveals a nodal d-wave 
type excitation spectrum inside the pseudogap, testifying to 
the shared origins with the superconducting state. Importantly 
however, the observed quantum diamagnetism [ 6 ] terminates 
at very low but finite underdoping and thereupon the magnetic 
signal is dominated by spin response. This is in line with the 
ARPES data on BSCCO El, infrared ellipsometry on RBCO 
ifTOl and the most recent STM data fTTl . which are suggestive 
of two pseudogaps: the nodal one, reminiscent of a supercon- 
ductor and decreasing at extreme underdoping, and the antin- 
odal one, which appears to be dominated by spin correlations 
and remains large. 



Two kinds of new results are reported: First, guided by 
broadly accepted microscopic features of cuprates, a theory 
for the long-distance charge 2e sector of the pseudogap state 
is constructed and shown to provide quantitative understand- 
ing of the observations in O, including the measured upper 
critical field at T = and the quantum vortex liquid and 
solid at lower fields. Furthermore, a dSC dome in the doping - 
magnetic field-temperature (x-H-T) phase diagram is found 
to be enveloped by a larger, charge 2e Cooper pairing dome, 
dominated by quantum fluctuations of the superconducting or- 
der parameter This larger dome collapses to T — at finite 
x = xq ~ 0.01, while the amplitude of the microscopic spin- 
singlet pairing term Ajk remains large as x — > 0. Between xo 
and x c ^ 0.055, the ground state is the quantum disordered 
dSC, the order in \I/ (~ (e %9jk )) preempted by free vortex- 
antivortex excitations of \Ajk\e L ° jk . Two distinct pairing en- 
ergy scales - for the charge and A for the spin sector - 
provide a natural explanation for "small" and "large" pseudo- 
gaps observed by various experimental probes |[Tl[9l [T0l[TT1l . 

Second, a sharp distinction is drawn between the dSC 
(off-diagonal or particle-particle) and particle-hole (diagonal) 
regimes of the microscopic theory in terms of the symmetry of 
the action for the bond phase 6jk of Aj&: in the former case 
this symmetry is of the XY type, while in the latter it is re- 
lated to a compact gauge symmetry, and thus much larger. A 
quantum phase transition between the two regimes is likely at 
extreme underdoping x < xq. The XY regime is dominated 
by the center-of-mass motion of the Cooper pairs, while the 
relative motion fuses into low energy physics in the "compact 
gauge" regime, with fundamental consequences. By analogy 
with its more familiar s-wave (bosonic) cousin l4l 171 [T2l [T3lL 
we dub the above sequence of events - unleashed by inten- 
sifying quantum fluctuations of Ojk on approach to the Mott 
insulating state - a "d-wave duality". 

The starting point is a general strongly correlated micro- 
scopic Hamiltonian which almost certainly underlies the es- 
sential physics of cuprate superconductors: 

& = - E 4*v + E ^MAi - 4A\ + ( h - c -) 

+ -j- E AlAj + E Up^nn + E Wifiinj ■ (!) 

In c\ a (c icr ) creates (annihilates) electrons of spin a on site 
i of the Q1O2 lattice, hi G = c\ a Ci a , hi = + 77^, i* is an 
effective hopping, A^ is the Hubbard- Stratonovich "pairing" 
operator, the integration over which yields the superexchange 
term J e ^ J2(ij) [&i ' an d U p is an arbitrarily large 



2 



repulsion, a purely mathematical tool introduced to suppress 
double occupancy iTH. {W^ } are the extended-range (near- 
est, next-nearest,...) Coulomb and possibly phonon-induced 
interactions. At half-filling (x = 0), hi = 1 and all but 
the "pairing" terms in |TJ) drop out; the result is the S = \ 
Heisenberg Hamiltonian and the ground state is the Neel-Mott 
antiferromagnet (AF). At x > 0, the precise nature of the 
ground state of strongly interacting Hamiltonians ([I]) is un- 
known but we assume a domain of its parameters where a 
correlated BCS-type d-wave superconductivity is stabilized, 
as established in a myriad of experiments Q. 

We anchor our analysis to the off-diagonal saddle point of 
{!]): \Ajk\ = Ad = ±A at x(y) bonds. As x decreases toward 
x c and beyond, the increasing quantum disorder in 6jk pro- 
gressively destabilizes the dSC ground state. Nevertheless, we 
persist in following this progression all the way to the Heisen- 
berg AF at x = 0, which, in a sense precisely defined below, 
is an "infinitely" strongly quantum fluctuating dSC. 

To follow this path of growing fluctuations in 6jk, we 
integrate out the fermions in ([T]) about the above saddle 
point. This is (much!) easier said than done: the pres- 
ence of gapless nodal fermions in a d-wave superconduc- 
tor makes the resulting action severely non-analytic. Still, 
as long as we stay focused on physics intrinsically tied to 
Cooper pairs and steer clear of nodal quasiparticles (or gap- 
less spinons discussed in 0), the following (bond) phase- 
only action for the correlated dSC captures the essentials: 
J V(fi exp(— J drC d XY ), where 



qd 

&XY 



- Jl ^ COs(<£i - <pj) ~ J 2 ^ COS (^i - Vj) 
rnnn bnnn 

- COs(D#) + £nodal[cOs(<Pi - (fj)} + £ C ore • (2) 



In ([2]), {(fi} are identical to the bond phases {Ojk} of {Aj^}, 
with subscript i referring to the center of each bond (jk), the 
first term is the charge 2e Berry phase, nn and r(b)nn are 
nearest and (inequivalent!) next-nearest neighbors on the {i} 
lattice, □ indicates plaquettes of the CuC>2 lattice, and DO = 
#12 — #23 + #34 — #4i around each plaquette. The meaning of 
([2]) and its relation to ([I]) are detailed in [ 7 ] and in Fig. 1. We 
will use ([2]), and lean on ([T]) when necessary, to analyze the 
low-energy physics of cuprates. 

To extract the long distance effective theory of quan- 
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FIG. 1 : A sampler of processes in |TJ contributing to various cou- 
plings in ([5]). The vector potential A em couples only to hopping t* . 



Upon introduction of the complex field = II ^ 
pled to e~ l(fii , the integration over {(pi(r)} yields 



iQi cou- 



dr 



D* i (r)cxp[/ ^-^Bm(T)\ 2 )- 

i 

r^r r ^^Kr^TlQoiKujUTf^jiT')} (4) 

J 'V J 'V 



as the representation of ^ in terms of {\E^(t)}. In 
B(\z\ 2 ) = \n[J^dcpexp(ze-^ + z*e^)] and Q^^oo) 
is a 2 x 2 matrix, with diagonal elements J\ cos(k x a) + 
J2 cos(/c y a) + ^o^ 2 and J2 cos(k x a)-\- J\ cos(k y a)-\-^tvoiJ 2 , 
and off-diagonal ones J[cos((k x — k y )a/2) + cos((/c x + 
ky)a/2)]. The 2x2 form of Q^fcu) simply reflects 
Ji 7^ J2 and the resulting inequivalence of the nearest neigh- 
bor sites of the {%} lattice, centered on bonds of the original 
Cu0 2 lattice (Fig. 1). 

Several simplifications are built into ^ in anticipation of 
our next step: K ^ is assumed to be small in the regime 
of interest (see later text, however). Furthermore, at long 
distances relevant to 0, we can neglect the Berry phase in 
([2]) which induces 2e superstructure modulations in CuC>2 lat- 
tice at intermediate lengthscales. r v , the vortex tunneling 
time, is a short time cutoff: r v ~ m v a 2 , where m v is the 
(anti) vortex mass, equal to only a few electron masses in un- 
derdoped cuprates 171 [T5l. r v should not change dramatically 
within the window of x probed in [6]. 

Assuming |\E^(t)| 2 <C 1, focusing on long distances and 
low energies, and restoring H via V — > V — z(2e)A gives 
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/ d(T/T v )d 2 (r/a)C d QGL [V(r, r)], where 
*| 2 +7|(V-i(2e)A)*| 2 +a|*| 2 - 



*| 4 , (5) 



turn fluctuations in {e l 
Hubbard- Stratonovich fields, {11^ (r)} and {Q 
pie terms like exp( J -\ cosf^ — tPn)) in 



i(r)h 

{iJ} cos(ifi - ipj)) inQas 
i cos ifi + Qi sin 

i ^ i 

X)(n<,Qi)(i|(J cosfeaaJ-^jXn^Q 



7r = K /8Jjrl 7 = a 2 J/SJ 2 r 2 , a = 
e^ k }, we introduce two r^a/ J + \ J\ + ^ J 2 , and J = J + Ji + J; 



to decou- 



(3) 



a=x,y 



(1/2 J d r v ) - 1, J d = 
. Note that, consistent 
with our strategy, £ no dai @ produces only higher (third order, 
non-analytic) derivatives in ([5 ). 

The physical meaning of (5b is important: Cq GL repre- 
sents the low energy, long-distance charge 2e sector of ([T]) 
in the regime of its parameters and doping characterized by 
strong quantum superconducting fluctuations. While the pro- 
cedure that led us to Q from ([2]), via ([3] [4]), bears a superfi- 
cial resemblance to the standard Gor'kov's derivation of the 
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Ginzburg-Landau functional from the BCS theory, the actual 
physics is different: by setting ^ to a complex constant in the 
standard Gor'kov- Ginzburg-Landau functional and focusing 
on the quadratic term, one is probing the weak-coupling in- 
stability of the electronic system toward formation of Cooper 
pairs, i.e. a finite complex A, accompanied by full phase co- 
herence. Since the pairing susceptibility of a typical normal 
system diverges as T — > 0, such instability is always present 
as long as there is an attractive coupling constant in a given an- 
gular momentum channel. In contrast, ^ in Q is testing for 
phase order in e t9jk in the system in which large | Ajk | 7^ is 
already established. The susceptibility to such phase order is 
typically finite, something the reader can easily check by us- 
ing ([3]) and applying our procedure to the simple quantum XY 
model lfT6l . This fact of much theoretical importance consti- 
tutes the foundation of various dual descriptions of strongly 
correlated systems (7] [12] [13|: a single effective action ^ 
describes both ordered and quantum disordered ground states 
- depending on the strength of its coupling constants - just 
like the original microscopic Hamiltonian. Obviously, ([5]) as 
it follows from ^ is only a coarse-grained caricature of ([I]). 
Nevertheless, its key features faithfully reflect general proper- 
ties of the underlying microscopic theory. 

The first among these are Mott correlations which, as x — > 
0, suppress double occupancy and drive kinetic energy toward 
zero. This implies Jd(x) — > 0: it becomes impossible to es- 
tablish phase coherence between on different bonds (Fig. 
1). This is reflected in a(x — > 0) — > +00 §5§. In the lan- 
guage of ([5]) x = acts as an infinite "temperature" and thus 
the familiar Heisenberg AF at half-filling is "infinitely far" re- 
moved from a dSC (and vice versa!). This statement is made 
more precise from the standpoint of the microscopic theory 
(p} later in the text. 

While the exact calculation of t*(x), Ad(x), and Jd(x) is 
beyond reach, the above trend is manifest in various approx- 
imations (t* ~ xt, A d - J eff as x -> E [13 HI). In 
any case, the precise form of Jd(x) is not needed for our 
considerations: we can expand a(x) w c*o(xo — %) about 
a(x = xq) = (l/2Jd(xo)r v (xo)) — 1 = 0, the above argu- 
ments clearly implying that xo is finite. A reasonable estimate 
gives xo ~ Ad/t, where t is the bare hopping. 

With these simplifications, we use ([5]) to address certain ba- 
sic features of the cuprates' x-H-T phase diagram. In uniform 
H , we follow Abrikosov and expand the quadratic terms in ([5]) 
in charge 2e Landau levels (LLs) (distinct from single electron 
LLs!) tf (r, r) = £\ ^ (r, r), j = 0, 1, 2, . . . : 

CO -j 

£gl - £ W*il 2 + a^H^tf] + , (6) 

3=0 

where aj = ao (xo — x + (2j + 1)H/H' c2 ) . aj=o(#, H) = 
determines the upper critical field 

H c2 (x,T^ 0) = ^ 2 (x - x ) , H f c2 = a /47re 7 . (7) 

In Fig. 2(a) we compare ^ to (6) and find good agreement 
with the choice xq « 0.01 and if£ 2 ~ 9.4 Tesla/% of doping. 



The upper critical field ([7]) is a mean-field concept - it 
decribes the line in the H-x phase diagram along which 
<Sg GL [\I/(r, t)) (pi) has an absolute minimum for a single 
co nfiguration o f w(r,r) = ^ m f(r). At H = 0, \P mf = 
y^2ao(x — xo) 7^ 0, implying {e l9ij ) 7^ 0, for x > xq. For 
x < xo, however, the minimum of Sq GL is at \£(r, r) = 
{e l0i o) = 0! Thus, for low x, the above mean-field theory of 
Sxy Pf and, by inference, of its microscopic parent Hamilto- 
nian is intrinsically different from Z?6tf/z the standard BCS 
theory and the slave boson mean-field approximation and the 
related Gutzwiller-oroiected BCS wavefunction |[3l [T7l[T8l . 




FIG. 2: (a) H [Tesla] vs. x phase diagram of LSCO adapted 
from |6|. All data points are from |6|. Thick (blue) dashed line is 
H c2 (x,T -> 0) |7j, with x = 0.01 and H' c2 ~ 9.4 Tesla/%. 
H m (x,T — > 0) is the boundary of the dark (red) region. H m (x) 
turns convex by T ^ 4K, in line with ( |11|12| ). The lighter (blue) 
region shows increasing spin response |6|; (b) x-T phase diagram of 
cuprates based on \2\5\ . Note that T m f is roughly of the order of the 
"small" charge 2e pseudogap ^ m f Ad while Ta ~ T* corresponds 
to the "large" pseudogap Ad. 

This difference stems from an important physical informa- 
tion captured by ([5]): In a BCS-style weak-coupling theory, it 
is the kinetic energy part of ([I]) that is large at the point of the 
quantum phase transition, while the amplitude A^ — > 0. Thus, 
the moment the self-consistent solution A^ 7^ appears, the 
phase coherence between different bonds is automatically es- 
tablished. \I/ m f is immediately finite and proportional to A^. 
Incidentally, this is the likely physics behind the closing of the 
dSC dome on the overdoped side of the phase diagram in Fig. 
2(b), at x = x wc . Beyond x wc , there are no superconducting 
fluctuations in the ground state; the system is in its underlying 
normal state, most likely the Fermi liquid. 

In the underdoped regime, near xq, the situation is reversed. 
At low x, we are dealing with rapidly growing Mott correla- 
tions and the ensuing suppression of coherent motion of 2e 
charge. Now, it is the kinetic energy that is becoming small 
while the spin-singlet pairing |A| remains large. However, the 
superconductivity now does not appear the moment t* 7^ - 
in a strongly correlated system the kinetic energy must rise 
to a finite fraction of A^ before even a mean-field solution 
^mf 7^ becomes possible at x = xo. This crucial fea- 
ture of the microscopic theory ([T]) is faithfully captured by 
Eq. ([5]) and the mean-field approximation \£(r, r) — > \I/ m f. 
For x < xo, even though A^ remains large (~ Jeft), ^mf 
collapses to zero (Fig. 2), reflecting enormous quantum flue- 
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tuations in {%}, unrestrained by the coherence between dif- 
ferent bonds. In contrast, the slave boson and the Gutzwiller- 
projected mean-field theories [17] [T8l produce a dSC the 
moment kinetic energy, ~ xt, becomes finite. 

How are the above mean-field arguments generalized to 
T / 0? Various coefficients in ^ have temperature depen- 
dences that are difficult to determine directly from a strongly 
interacting theory ([I]). Again, a simple Gaussian correction to 
the quadratic term in ^ aj — > aj (x, H, T) w a (#o — # + 
(2j + 1)H/H' c2 -f T 2 /0 2 ) will suffice for our needs and gives 



i7 c2 (x,T)=^ 2 (x-xo-T 2 /e 2 ). 



(8) 



A corollary of (|8j) is that, at H = 0, T mf = 6^ - < 
Ta ~ |A| defines the line in the x-T phase diagram below 
which ^ m f(x, T) = ^/2ao (x — xq — T 2 /6 2 ) becomes finite 
(Fig. 2(b)). The estimate from the underdoped side of "Nernst 
dome" yields 6 - 37K/V% in LSCO. 

The above mean-field transition to \P m f(r) / from 
^mf(r) = is only a crossover: Eqs. ( |7|8| ) define the un- 
derdoped side of the "Nernst dome" (or "diamagnetic dome") 
in the x-H-T phase diagram [ 6] along which enhanced dia- 
magnetism and (anti)vortex fluctuations become observable 
(see Fig. 2). These are but a couple of manifestations of the 
rise of superconducting correlations and of |\£(r, r)| in ^ 
sampling large local values even in absence of the long-range 
off-diagonal order. Actually, much is known about the fluc- 
tuations in ([5] [6j (see [19] |20) and references therein). At 
H = 0, the quantum phase transition to the true supercon- 
ducting ground state is shifted to x c from xq\ x c — xq ~ Gi, 
where Gi is the dimensionless coefficient of the quartic term 
in ^ |19]; Gi - 0.045 from our fit to 0. At finite T, the 
true superconducting transition line emanates from x c as 



T c (x) 



a: 



z/(z-2) 
Xd 



~ \X - 



, (9) 



where p s (x), X s (x), Xd(%), £sc(#)> v -> and z are me superfluid 
density, penetration depth, diamagnetic susceptibility, corre- 
lation length, its exponent, and dynamical exponent, respec- 
tively, all associated with the quantum superconducting transi- 
tion at x c . Eq. ^ predicts v = vxy ~ 0.667 and z = 1 but a 
more detailed treatment |7], incorporating the Berry phase ^ 
can lead to different values. The regime between the mean- 
field crossover and the true transition is dominated by quan- 
tum vortex-antivortex fluctuations (Fig. 2(b)) |2T1 . 

Similarly, the true transition for H ^ in §5§ takes place 
only far below H c2 (x,T) ([§}, at H m (x,T). The x-H-T 
phase diagram between H C 2(x,T) and H m (x,T) is occu- 
pied by the quantum vortex liquid (Fig. 2(a)), exhibiting 
strong vortex fluctuations and enhanced diamagnetism. Along 
H = H m (x,T), this quantum vortex liquid freezes into a 
solid, and any weak pinning ensures superconductivity. The 
shape of the H = H m (x, T) line depends on whether one is 
in the "low" or "high" field regime of ^ [19]. For H > H b , 
the quartic term in ([6]) is ineffective in mixing the LLs of the 
quadratic part; in the opposite limit, H <C H b , this mixing 
is strong and the H ^ behavior is dictated by the quantum 



vortex-antivortex unbinding at H = [19J |20l. An estimate 
OH gives H b ~ (Gi/13)i^ 2 - 4 Tesla in LSCO, about mid- 
way through the fluctuation region explored in (6). 

At high fields, H ^> H b , ^ is dominated by the low- 
est Landau level (LLL) (J = 0). Sq GL 
Jd(r/r v )d 2 (r/a)C d QGL _ LLL [^o(r,r)},wh^: 



°QGL-LLL 



£gl-lll =7r|^o| 2 +^o(^^T)|^ | 2 + ^|^o| 4 , (10) 

with a j=0 = a (x - x + H/H' c2 + T 2 /6 2 ). The func- 
tional integral J V^/ exp(— Sq L _ lll ) is confined to the LLL, 
* e tfo(r,r) = ®(T)Ui( z ~ ^(r))exp(-|z| 2 /4), where 
z = (x + iy)/£, £ = ^l/2eH and {zi(r)} are the quan- 
tum zeroes (vortices) of ^(r, r) fT9ll . A prominent feature 



of ( 10 ) is that the physical quantities computed from it exhibit 
single-parameter scaling: ground state energy, magnetization, 
low T thermodynamics, etc., are all universal functions of 
g(x,H,T) = const, x (x - x + H/H f c2 + T 2 /@ 2 )/H 2 ^ 
only. Various scaling functions are computed in fT9l . In par- 
ticular, the quantum vortex liquid-solid transition in {^(r)} 
observed in [ 6 ] takes place at g(x, H, T) = g m < 0, where 
g m is a universal number, g m « —3. This leads to 



H m (x,T)~(x-x -T 2 /G 2 ) 3 / 2 . 



(11) 



H C 2(x, T) (|7j) corresponds to g = 0; thus, the fluctuation re- 
gion in this regime expands to H C 2 — H m ~ Gi 2//3 . 

In the opposite limit, H <C H b , the H = quantum tran- 
sition governs scaling properties: the universal scaling func- 
tions depend on (, sc (x, T = 0)/£ and H m (x) emerges from x c 
as 



H m (x) ~ (x 



(12) 



Note that both ( fTT| ) and $12) are convex lines while the ob- 
served H m (x, T — > 0) is concave (Fig 2(a)). The likely cause 
is pinning of quantum vortices by the Cu02 lattice (or disor- 
der), absent from continuum theory ([5]). This explanation is 
supported by the H m (x) curves in [ 6 ] turning convex at a fi- 
nite, but still low, T (~ 4 K <C |A|) suggestive of thermal 
depinning and the H m (x) lines for different T all terminating 
at the same x c « 0.055 as H — > 0. The number of data points 
in (6) is insufficient for a detailed fit to ( |11|12| ) - it would be 
interesting to test these predictions in further experiments. 

The physics captured by §5§ and discussed above reflects 
the center-of-mass motion of Cooper pairs. To display this 
more clearly consider slowly changing phase Ojk, far from the 
cores of vortex or anti vortex defects. Since electrons live on 
the sites of Cu02 lattice it is convenient to define two sets 
of site variables: phases {<j)j} as e 1 ^ = J2k e%9jk /\ J2k e% ° jh I 
where k runs over four bonds emanating from site j, and com- 
plex numbers {s w } as s w = e i012 - e i023 + e i034 - e^ 41 
where uo denotes a "virtual" site at the center of each (1234) 
plaquette of the Cu02 lattice. Up to and including sec- 
ond order derivatives one can invert these definitions to find 



i if) 

where e o k 
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l22l and O v - k [{0j}, {scj}] is the "irreducible" part of Ojk which 
cannot be reduced to site phases. By inserting the above ex- 
pression for Ojk into {!]) and simply repeating the steps that led 
from ^ to (5]), one arrives at precisely the same phase struc- 
ture as in ([5 ), except that now 0(r, r) and V0(r, r) — (2e) A 
replace the corresponding derivatives of Evidently, (j)j 
plays the role of the overall center-of-mass phase of a Cooper 
pair. Importantly, 6 r - k does not appear at the second order level 
in derivatives and is naturally interpreted as corresponding to 
the relative motion of the paired spin-singlet bonds in ([I]). 

This brings us to the second feature that transpires from §5§ 
and concerns extreme underdoping x < xq: in this regime 
a — > +00 and off-diagonal correlations are maximally sup- 
pressed. Such limit of extreme quantum fluctuations is dif- 
ferent for the d-wave bond phase Ojk than for the site phase 
in an s-wave (bosonic) XY duality (I2j [T3ll . In the s-wave 
case, the dual opposite of a superconductor is a Wigner crys- 
tal of Copper pairs. This remains an entirely center-of-mass 
affair, with the quantum disorder in the corresponding site 
phase promoted by suppressed density fluctuations of posi- 
tionally ordered Cooper pairs. But off-diagonal correlations 
are still strong, as exemplified by quantum diamagnetism and 
vortex-antivortex fluctuations, and the ground state is still well 
described by the XY-like phase action; these properties are 
broadly similar to the intermediate "Quantum vorticity" state 
in Fig. 2(b) although there are important differences as well, 
like the presence of spinful nodal fermions in the latter. 

By contrast, in the d-wave case, the extreme quantum dis- 
order limit (x — > 0) of the action ^ does not have the XY 
symmetry: t* — > and the XY parts of ^ disappear, leaving 
behind only the plaquette-type action, since terms like K in- 
clude only A's (Fig. 1). This demise of the XY behavior in 
([2]) and ([I]) - occasioned by Mott suppression of kinetic en- 
ergy - translates into collapse of off-diagonal correlations and 
is identified here as the microscopic physics behind the ob- 
served "two pseudogaps" in underdoped cuprates l9l [T0l[T71l . 
as alluded to earlier in the text (Fig. 2(b)). To see this explic- 
itly, we consider below the off-diagonal, "anomalous" elec- 
tron propagator following from {!]): F(r — r',r — r') = 
(cj(r, r)c|(r / , r')). Finite F implies superconducting off- 
diagonal long-range order. 

We can actually compute F with the accuracy needed to 
make contact with the recent experiments (U [lOl [TTJ. First, 
within the "Quantum vorticity" state in Fig. 2(b) it suffices to 
keep only the leading relevant derivatives and we can replace 
e i0 jk _^ e i0™ = ( e i0j + e <0fc)/|e^ + e i<t>k \, as explained 
above. Next, we define new fermions f icr as c icr = fi G e^ i , 
cj = e~^ i f\ . This leads to: 



r ,r 



T')K(*)(f}(r,T)fl(r\r')) , 



(13) 



accurate up to leading derivatives, where is nothing but 
the expectation value of the superconducting order param- 
eter computed from ([5]). The precise structure of the off- 
diagonal /-fermion propagator ( 13 ) is a complex problem, as 
/ fermions behave as neutral d-wave quasiparticles interact- 



ing via gauge fields generated by spacetime derivatives of the 
phase \(\> Q; nevertheless, for the purposes of overall ener- 
getics and at the level of resolution available in |5] [TO] [TTJ, 
such propagator behaves simply as - a somewhat smeared 
version of - the ordinary BCS d-wave anomalous propaga- 
tor. Furthermore, Refs. (9J [TOl [TT1 are sensitive only to short 
range off-diagonal order and thus should be replaced by 
\I/ m f . After the Fourier transform we finally obtain: 



(v F -(k-k F )) 2 



|A,| 2 ' 



(14) 



where A^. = Acos(2(/>k), </>k is the angle going around the 
Fermi surface and vp and kp are the Fermi velocity and 



momentum, respectively. From (14) it follows naturally that 
A should be interpreted as the "small" pseudogap, mea- 
suring charge 2e off-diagonal correlations and collapsing near 
xq, while the spin- singlet pairing gap A itself remains large 
(Fig. 2(b)). 

The above collapse of "small" pseudogap \I/ m f A at x < xq 
implies that vortex-antivortex pairs, the staple of XY-type 
physics, are displaced as relevant excitations. To see what 
replaces them, we transform c iy(T — > ad\ _ a on one of the 
sublattices of the CuOolattice, turning the spin-singlet pair- 



ing |A|e^ fe ct ?a 4 5 _ a M into hopping |A|e*"ct a d fc|<r . Af- 
ter this transformation, {±0jk} on alternating bonds are cast 
in the role of a gauge field coupled to the staggered 

"charge", switching from +1 to —1 between the two sublat- 
tices. In this language, the plaquette terms in ^ are nothing 
but the action for {djk}, invariant under the compact gauge 
transformations — > + Q — This is a much larger 
symmetry than the XY one of ^ at finite x and t* . The rel- 
evant excitations of a compact gauge theory are monopoles 
and antimonopoles in {cijk} - in turn, monopoles in 
are closely related to the relative motion of Cooper pairs and 
the bond phases {0j k } defined earlier. This is easily appreci- 
ated by observing that the plaquette term in ^ always equals 
±1 when Ojk is restricted to 0™ alone. Only by allow- 
ing unrestricted fluctuations of Ojk and including the relative 
phases v - k , can one successfully proliferate monopoles and 
antimonopoles in a^. With ([TJ serving as the microscopic 
model, the monopoles will be in their plasma phase and thus 
cofinining. Actually, despite apparent mathematical opaque- 
ness of our argument, this is just one among many disguises 
of the traditional Heisenberg-Mott AF Q. This opaqueness 
is an unavoidable consequence of our persistence in follow- 
ing the fate of a dSC all the way to the half-filling, the route 
along which the quantum fluctuations in the phase Ojk of its 
gap function progressively become enormous and ultimately 
unrestrained. In this sense, it is the standard Heisenberg-Mott 
AF that stands at the far end (x = 0) of the "d-wave dual road" 
for correlated dSC l23l . 

How is the XY symmetry recovered at finite x and t*7 The 
J\ term in ^ provides an illustration (Fig. 1): it moves charge 
2e "Cooper pair" from one spin- singlet bond to the next, a 
process vital for establishing local superconducting phase co- 
herence. In the gauge theory language this is a Higgs term: 
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cos(Oij — Oki) — > cos(dij + a,ki) \ Thus, the charge 2e motion 
and the associated kinetic energy correspond to the symmetry 
breaking Higgs terms in the gauge theory dialect. The pres- 
ence of Higgs terms suppresses free monopoles in {a^} and 
promotes its gauge sector to physical reality - the result is the 
conservation of vorticity in {%} and the XY symmetry of 
its action ([2]). Thus, Copper pairs "emerge" from the under- 
lying dynamics of ([T]) within the "Quantum vorticity" region 
in Fig. 2(b), as dual partner of well-defined vortex-antivortex 
excitations. One is tempted to speculate that a quantum transi- 
tion between the "compact" and "XY" regimes takes place at 
finite x ~ xo, with the AF+x ground state(s) in Fig. 2(b) reg- 
ulated by monopoles and antimonopoles in ajk and reflecting 
the relative motion, instead of vortices and antivortices asso- 
ciated with the Cooper pair center-of-mass fluctuations in 9jk- 
Such ground states are the true "competing orders" of a dSC, 
arising from the particle-hole (diagonal) channel, and likely 
include non-uniform AF and stripes [201. 
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